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Abstract 

Non-Gaussianities of the primordial density perturbations have emerged as a very powerful possible 
signal to test the dynamics that drove the period of inflation. While in general the most sensitive 
observable is the three-point function in this paper we show that there are technically natural 
inflationary models where the leading source of non-Gaussianity is the four-point function. Using 
the recently developed Effective Field Theory of Inflation, we are able to show that it is possible 
to impose an approximate parity symmetry and an approximate continuos shift symmetry on the 
inflaton fluctuations that allow, when the dispersion relation is of the form u ~ c s k, for a unique 
quartic operator, while approximately forbidding all the cubic ones. The resulting shape for the 
four-point function is unique. In the models where the dispersion relation is of the form oj ~ k 2 /M 
a similar construction can be carried out and additional shapes are possible. 



1 Introduction 



In the last few years, it has become more and more clear that single-clock inflation, i.e. 
inflationary models where there is only one degree of freedom driving inflation, can produce 
large and detectable non-Gaussianities. While standard slow roll inflation cannot produce 
large non-Gaussianities [1], several models in single field inflation have been proposed that 
result in a detectable level of non-Gaussianities [2, 3, 4, 5, 6, 7, 8]. Furthermore thanks 
to the development of an Effective Field Theory for the Inflationary perturbations [6, 7], it 
has become clear that the Lagrangian for the fluctuations can incorporate large interactions 
without spoiling the naturalness of the quasi de Sitter background. 

Non-Gaussianities represent a powerful signal about inflation for two reasons. First con- 
trary to what happens for the two-point function, non-Gaussianities are sensitive to the in- 
teracting part of the inflaton Lagrangian which is clearly more interesting than the free-field 
part. Second, the symmetries of the inflationary spacetime force the two-point function of the 
cosmo logical fluctuations to be quasi scale invariant. This leaves freedom only to vary its am- 
plitude and to add a small scale-dependence as we vary the inflationary models. This means 
that measurement of the inflationary two-point function reduces in practice to measuring only 
a couple of numbers. On the other hand, the same symmetries of the inflationary spacetime 
are not able to constrain as much correlation functions of higher order. The possible signals in 
this case are continuous functions of several variables, very similar to the angular dependence 
of a scattering amplitude [7, 9]. This is clearly a much more powerful signal. 

So far, the observational search for non-Gaussianity as originated from the period of infla- 
tion has concentrated on the three-point function. Several models predicted such a signal as 
the leading source of non-Gaussianity. The purpose of this paper is to show using the Effective 
Field Theory of Inflation [6] that there are technically natural models of inflation where the 
level of the (connected) four-point function is much larger than the three-point function, and 
where the leading source for detection is indeed in the four-point function. Using the Effective 
Field Theory of Inflation we show that it is possible to impose on the inflaton fluctuations 
both an approximate continuous shift symmetry and an approximate parity symmetry that, 
for models not close to the de Sitter limit, allow for a unique large quartic operator, 7r 4 , where 
7r represent the inflationary fluctuations. These symmetries approximately forbid all cubic 
terms. The resulting shape of the four-point function is unique. A similar construction is 
carried out for models near the de Sitter limit where more that one shape becomes possible. 

2 Effective Field Theory of Single-Clock Inflation 

In this section we briefly review the effective action for single-clock inflation. This effective 
action was developed in [6, 15] and we refer the reader to those papers for a detailed expla- 
nation. The construction of the effective theory is based on the following consideration. In 
a quasi de Sitter background with only one relevant degree of freedom, there is a privileged 
spatial slicing given by the physical clock which allows us to smoothly connect to a decelerated 
hot Big Bang evolution. The slicing is usually realized by a time evolving scalar (f>(t), but 
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this does not necessarily need to be the case. To describe perturbations around this solution 
one can choose a gauge where the privileged slicing coincides with surfaces of constant t, i.e. 
S(j)(x, t) = 0. In this 'unitary' gauge there are no explicit scalar perturbations but only metric 
fluctuations. As time diffeomorphisms have been fixed and are not a gauge symmetry any- 
more, the graviton now describes three degrees of freedom: the scalar perturbation has been 
eaten by the metric. One therefore can build the most generic effective action with operators 
that are functions of the metric fluctuations and that are invariant under the linearly-realized 
time-dependent spatial diffeomorphisms. As usual with effective field theories, this can be 
done in a low energy expansion in fluctuations of the fields and derivatives. We obtain the 
following Lagrangian [6, 15]: 



= Jd 4 x^j ^M^R + M^Hg 00 -M^(3H 2 + H) + 



+ ^M 2 (t) V° + I) 2 + ^M 3 (t) V° + I) 3 + 



where we denote by SK^ U the variation of the extrinsic curvature of constant time surfaces 
with respect to the unperturbed FRW: 5K^ U = K^ u — a 2 Hh^ u with h^ u being the induced 
spatial metric, and where M 2j 3 and M 1)2 ,3 represent some time-dependent mass scales. 

Let us comment briefly on (1). The first term is the Eistein-Hilbert term. The first 
three terms are the only ones that start linearly in the metric fluctuations. The coefficients 
have been carefully chosen to ensure that when combined the linear terms in the fluctuations 
cancel. The action must start quadratic in the fluctuations. The terms in the second line start 
quadratic in the fluctuations and have no derivatives. The terms in third line represent higher 
derivative terms. Dots represent operators that start at higher order in the perturbations or in 
derivatives. This is the most general action for single field inflation and in fact it is unique [6] . 

The unitary gauge Lagrangian describes three degrees of freedom: the two graviton he- 
licities and a scalar mode. This mode will become explicit after one performs a broken time 
diffeomorphism (Stiickelberg trick) to reintroduce the Goldstone boson which non-linearly 
realizes this symmetry. In analogy with the equivalence theorem for the longitudinal com- 
ponents of a massive gauge boson [19], the physics of the Goldstone decouples from the two 
graviton helicities at high enough energies, equivalently the mixing can be neglected. The 
detailed study of [6, 7] shows that in most situations of interest this is indeed the case and 
one can neglect the metric fluctuations 1 . 

As anticipated, we reintroduce the Goldstone boson (tt) by performing a broken time-diff., 
calling the parameter of the transformation — it, and then declaring n to be a field that under 
time diff.s of the form t — > t + £°(x) transforms as 

ir(x) ->■ n(x(x)) = n(x) - (j°(x) . (2) 

In this way diff. invariance is restored at all orders. For example the terms containing g 00 in 



1 Equivalently, the neglected effects are suppressed by slow-roll parameters or by powers of H/M P i. 
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the Lagrangian give rise to the following terms: 



9 °° - 9( L + / )3( LVV - 9 » + 2^+W (3) 

We refer to [6] for details about this procedure. If we are interested just in effects that are 
not dominated by the mixing with gravity, then we can neglect the metric perturbations and 
just keep the n fluctuations. In this regime, a term of the form g 00 in the unitary gauge 
Lagrangian becomes: 

^oo ^ _! _ 2tV - tt 2 + 4(^tt) 2 . (4) 

cr 

Further, we can assume that the n has an approximate continuous shift symmetry, which 
becomes exact in the limit in which the space time is exactly de Sitter [6]. This allows us to 
neglect to terms in n without a derivative that are generated by the time dependence of the 
coefficients in (l) 2 . Implementing the above procedure in the Lagrangian of (1), we obtain 
the rather simple result: 



-M 2 H ( n 2 - ^ + 2M$ (W 3 - *<*?n _ 4 M4 . a + 



' 2 \ " r " a 2 J 3 



, (5) 



where for simplicity we have neglected the terms originating from the extrinsic curvature as 
they are usually important only in a regime where the space time is very close to de-Sitter 
space [6]. In fact, notice that neglecting the terms in the last line in (1), we see that the 
spatial kinetic term is the form of HM^d^) 2 . In the limit H — >■ the spacetime approaches 
de Sitter and the coefficient of the spatial kinetic term goes to zero. If we consider the terms 
involving the extrinsic curvature in the last line of (1), upon reintroduction of 7r they give 
a spatial kinetic term of the form either (Mf + M 2 )(d 2 7r) 2 or HMf(diir) 2 . In general these 
terms are either higher derivative terms or suppressed by the Hubble scale, and so they are 
in general negligible at energies of order H. However, in the de Sitter limit the leading term 
HMp^diir) 2 goes to zero, and they become relevant. See [7] for a precise definition of when 
this is the case. We refer to this situation as the near de Sitter limit and we will explicitly 
describe the phenomenology in this limit at the end of the next section. 

We notice that when Mi is different from zero and we are not in the de Sitter limit the 
speed of sound of the fluctuations is different from one. We have the following relationship: 

Mi = . (6) 

In this limit there are two independent cubic self interactions, 7r(<9j7r) 2 and 7r 3 at this order 
in derivatives, which can induce detectable non-Gaussianities in the primordial density per- 
turbations. A small speed of sound (i.e. a large M 2 ) forces large self-interactions of the form 
7r(<9j7r) 2 , while the coefficient of the operator tt 3 is not fixed because it depends also on M 3 . 



2 Notice that this is not always the case. Interesting inflation models, both single field and multifield, 
have been recently proposed in which the tt fluctuations are protected only by an approximate discrete shift 
symmetry. See for example [10, 11, 12, 13, 8]. 



4 



By analysis of the cosmological data, one can therefore constrain (or measure) the parameters 
of the above Lagrangian. This approach has been recently applied to the WMAP data in [7] , 
giving constraints on M 2 and M3, as well as on the higher derivative operators that we have 
omitted in (5). This is the exact analogous of what happens for data from particle accelerators 
when the Precision Electroweak Tests of the Standard Model are carried out [17, 18]. 



3 A detectable four-point function from single field in- 
flation 

It is by now well established that single field inflation can produce a large and detectable 
level of non-Gaussianity in the three-point function. It is worth asking if it is also possible to 
have a large and detectable four-point function without at the same time having a detectable 
three-point function. The effective Lagrangian of single field inflation of sec. 2 is the ideal 
general set up to address this kind of questions. 

Restricting ourselves to the case where the Goldstone boson is protected by an approximate 
continuous shift symmetry, Ref.s [6, 7] show that in single field inflation there are only two 
ways to have a large three-point function: either by having a very small speed of sound c s 
for the fluctuations, or by the unperturbed solution being so close to de Sitter space that the 
dispertion relation of the Goldstone boson is of the form u 2 ~ k 4 /M 2 , where M is some mass 
scale related to M 2 ,3 and M 2 . 

In the case of a small speed of sound and away from the de Sitter limit, the large three- 
point function is induced by the operators 7r(<9j7r) 2 and fc 3 that are associated with the unitary 
gauge operators (Sg 00 ) 2 and (Sg 00 ) 3 . In particular, by estimating loop corrections reference [7] 
showed that if the coefficient of the operator (Sg 00 ) 2 is M| ~ HM^/c 2 using the c s < 1 limit 
of (6) then the operator (Sg 00 ) 3 is naturally of order M| ~ M^/c 2 (and viceversa). Both 
operators generate comparable three-point functions given by: 



(C 3 ) L 



(C 2 ) 3/2 L 2 



c 2 



V (7) 



where to estimate the effect we have taken the ratio of the cubic and the quadratic Lagrangian 
at energies of order H . It is customary to define 

<C^ 3 > = (2-) 3 * (3) (£i + h + k 3 )(2n) 3 (-^ W?) l 6 (8) 

— * — * 

where we have taken the equilateral limit of the 3-point function in Fourier space \k\\ = |/c 2 | = 
I £3 1 = k and 

(^ 1 ^) = (2 7 r 3 )^)(^ 1 + fc 2 )§, (9) 



so that 



InlC ~ 7^ 



=► /nl~^. (10) 
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The Planck satellite is expected to reach a limit on f NL ~ few [20] so we will limit 
ourselves to the case c s « 1. The reason why the natural value of M3 is M$/c 2 s can be 
quickly understood by doing the following manipulations of the single field Lagrangian 



d 4 x 



M 2 ,Hf. 2 2 {d t n) 



N&H Ad^f 2c 3 2 ■ 3 



(ii) 



where we have used (6), and we have redefined M| = c^M^/c 2 . We can perform a transfor- 
mation of the spatial coordinates: 



Ob ^ Ob X I Cg ■ 

and canonically normalize the field n 
to obtain 



dt d 3 x \f—g 




(12) 



(13) 



(14) 



(S\H\M^ 



1/2 



7T r 



c 3 



(S\H\M^ 



—7V 3 

1/2 "c 



where <9j = d/dx l . Notice that since we did not rescale time so the cutoff in energy can 
directly read off from (14). Since the kinetic part of the Lagrangian (14) is now Lorentz 
invariant, it is easy to read off the unitarity bound 



A 4 



<*\H\M&~ 



(15) 



In fact with this particular redefinition of the coefficient M3 in terms of £3 and c 2 we can see 
that for c 3 of order one the two operators n 3 and 7r c (<9j7r c ) 2 are suppressed by the same scale 
A 2 . Indeed it is straightforward to verify that loop corrections will generate the two operators 
with their relative coefficient c 3 of order one. We further notice that having a detectable 
three-point function means lowering the cutoff of the theory well below (iJMpJ 1 / 4 . Note that 
in terms of the cut-off the size of the non-Gaussianity is 



(C 3 ) 

(£2)3/2 



(16) 



as in fact they are generated by dimension six operators. 

Let us now consider the four-point function. The unitary-gauge operators (5g°°) 2 and 
(Sg 00 ) 3 contain many quartic operators. For example: 



M*(d t n) 



(17) 
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This operator induces a four-point function of the form 



(C 4 ) C 



r 4 



1 C 2 • (18) 



(C 2 ) 2 c 2 

The four-point function is usually parametrized by r NL = (C 4 )/(C 2 ) 3 so that, 

.2 ^4 



1 C 2 r NL ~ 1 . (19) 



c 4 c 4 



This value of tnl is tiny. The observational constraints errors on f^L and tnl scale as [22] 

A[W] - ~L , A[^lC 2 ] - , (20) 

pix pix 

where 7Vp ix represents the number of data points of the experiment. In order for a four-point 
function to be detectable, the value of r NL has to be a factor of about 10 5 larger than the 
value of Jnl allowed by the data [22]. Current limits from the WMAP satellite [7, 16] set 
Jnl ^ 10 2 , which from (10) implies that c 2 > 1CT 2 . Current limits on Tnl are therefore 
expected to be of the order 10 7 , and clearly the value in (19) is a factor of 10 3 too small. 
Even a detection of J'nl at its current upper bound of order 10 2 results in tnl ~ 10 4 — 10 5 
which will not be detectable by the Planck satellite. Planck is expected to produce a limit on 
tnl of order 10 6 . In the future 21-cm line experiments that could map a large fraction of our 
Hubble volume might reach f NL ~ 10" 2 , Tnl ~ 10 3 [23]. In order to detect the four-point 
function coming from these operators we will have to wait for quite some time and it could 
only happen if there is a detection of a large three point function 3 . So far we have considered 
only the quartic operator (<9j7r) 4 induced by the unitary-gauge operators (Sg 00 ) 2 and (5g 00 ) 3 . It 
is straightforward to see that the situation is the same also for the other operators associated 
with them. 

Notice that because of the non-linear realization of time-diffs, the coefficients of the quartic 
operators induced by (Sg 00 ) 2 and (5g 00 ) 3 where tied to the ones of the cubic operators. This 
was the reason why the induced four-point function was so small. However, one could imagine 
introducing the operator Mf(5g 00 ) 4 which starts quartic in the fluctuations and to arbitrarily 
set to zero the coefficients of the operators {Sg 00 ) 2 and (<5g 00 ) 3 .The fact that the coefficient 
of 7r is unrelated to the coefficients of the cubic terms had been already noticed in [24] 
for a subclass of the models we consider consisting of scalar field Lagrangians of the form 
P(((90) 2 ,0). However, without the identification of a symmetry protecting the generation 
of cubic terms, it is hard to imagine why one should have concentrated on the particular 

3 If one looks carefully at the current best limits on the three-point function presented in [7], one can notice 
that there is a small but non negligible region of the parameter space where values of c s as low as about 
10~ 2 are allowed. This is due to a partial cancellation of the three-point function in that region of parameter 
space, which however is not expected to happen at the level of the four-point function. This therefore implies 
a value of tnl of order 10 8 with such an high level of non-Gaussianity that could be detectable already by the 
WMAP satellite. Obviously, for this to happen we would have to be lucky enough so that inflation occured 
in that particular and small region of parameter space where the three-point function is suppressed. 
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Lagrangian allowing for a large quartic operator and small cubic ones. The point of our 
paper is that we are able to identify such a symmetry. 

Reintroduction of the field n would give rise to terms of the form 

M 4 W ) 4 M 4 4 (l67V 4 - 327r 3 (^7r) 2 + 247r 2 ( ( 9 M 7r) 4 -8^(^7r) 6 + ( ( 9 M 7r) 8 ) , (21) 

The term in 7r 4 induces a four-point function whose size is given by: 



(C 4 ) , 2 u 



Ml M 4 



C t-nl--^- . (22) 



Let us concentrate on the limit tnl ^> 1 (HMp X ^> A^) as this is the only regime where the 
signal is detectable. Upon canonical normalization, one can read that the unitarity bound of 
the theory is 

„ (jjgig (HMD ,„„ 

Av ~^r~^r- (23) 

tnl 1 translates into the condition HMp X ^> Afj, and we can rewrite tnl( 2 as 

.2 £4 



# 4 , X 

~ • ( 24 ) 



Thus by making the cutoff closer and closer to the Hubble scale we make the non-Gaussianities 
detectable. However in taking this limit we should worry about the fact that we have arbitrar- 
ily set to zero the operators (Sg 00 ) 2 and (Sg 00 ) 3 which we expect to be generated at loop level. 
Let us see at what level they are generated. Let us consider two of the terms in tt included in 
the operators (Sg 00 ) A : n 4 and 7r 3 (<9 M 7r) 2 . Upon canonical normalization, and neglecting factors 
of order one, we can see that these two operators appear in the Lagrangian as 



1 , 4 1 



-7T. 



At ' A6 T V2 



K 3 c (d,n c ) 2 . (25) 



U NL 

For t nl > 1, the operator n 3 , (<9 M 7r c ) 2 is effectively suppressed by an higher scale than 7r 4 . The 
same holds for the higher dimension operators that are included in (Sg 00 ) 4 : the higher the 
power in n they have the higher is the scale by which they are suppressed. This theory has 
an approximate and accidental ir — > —71 symmetry. 

If we now ask at what level these operators will generate the operators (Sg 00 ) 2 and (Sg 00 ) 3 
we realize that upon reinsertion of it these operator do not respect the symmetry n — > —n: 
they contain even and odd powers of ir. In fact their leading interaction term is odd. To 
generate them we will have to insert the operator 7r 3 (9 M 7r c ) 2 , effectively suppressing the level 
at which the operators (5g 00 ) 2 and (5g 00 ) 3 are generated. In other words, in the limit t nl > 1 
there is an approximate Z2 symmetry 71 — > —71 which is softly broken and protects the 
renormalization of the operators (Sg 00 ) 2 and (5g 00 ) 3 . 

For example, it is easy to estimate that the loop in fig. 1 induces an operator ir 3 suppressed 
by the scale (HM^) 3 / 2 . Since the structure of the Lagrangian is fixed by symmetries, by 
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diir 

7T 



Figure 1: Generation of a 3-point function from a loop diagram 



returning to the original normalization of the fields we conclude that the operators (Sg 00 ) 2 
and (Sg 00 ) 3 are generated with a coefficient of order HMp } : 

M*(6g°y HM^Sg 00 ) 2 ,HM^6g°y . (26) 

This is a very interesting result. In fact doing the same estimates that lead to (10) we 
conclude that the induced value of f^L is just of order one, while at the same time the t nl 
in (22) can be much larger than 10 5 . Therefore, we conclude that in single field inflation 
we can have a four-point function that induces a level of non-Gaussianities much larger than 
the three-point function and therefore can be detected even in the absence of a detection of 
a three-point function. Remarkably, this large four-point function is generated only by one 
operator: 7r 4 , and therefore there is a unique shape for it. Because it comes from derivative 
interactions this shape will be peaked on configurations in momentum space where all the 
momenta have comparable wavelengths. 

In particular, tnl can be larger than 10 6 , which is the approximate threshold for WMAP. 
So this signal could be waiting to be discovered in existing WMAP data. The theory and 
analysis of a four-point function in the CMB data is a largely unexplored area. We find that 
this result strongly motivates such a study. We will undertake it in a separate paper [21]. 

Of course, in order to have a large four-point function, it is not necessary that the operators 
(5g 00 ) 2 and (5g 00 ) 3 be exactly set to zero. It is straightforward to see that in the case of speed 
of sound different from one the condition for having the four-point function be the leading 
non-Gaussianity reads 

(#Mf.W2 

Finally, in models with a background space-time very close to de Sitter space [6, 7] the 
operator 7r 4 still gives a four-point function larger that the three-point function if the coef- 
ficient M 4 is large enough. There are however different possible shapes for the four-point 
function depending on the dispersion relation of modes with frequencies close to H . In the 
case where this is of the form oj = 'y 2 k 2 /M 2 , with 7 2 = (Mf + M 2 )/(2M 2 ), the operators of 
the form (SK^^) 2 and SK V ^SK** v dominate at horizon crossing, and the condition for having 
a four-point function larger than the three-point function reads 



K > . (27) 



M 25/4 M 3 M 13/4 

M 4 > 2 M A > 12 (oa\ 
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where for simplicity we have taken 7 < 1 in order not to have an effective speed of sound 
larger than one at energies below the cutoff. This estimate follows from realizing that the 
wavefunction in Fourier space for the canonically normalized ir c goes as 

1 M 2 1/2 . . 

TTc.fe ~ — F= ~ r~ , (29) 

which implies that at energies of order E, the size of the vacuum fluctuations of the operator 
ir(x) goes as 

1 / r m 2 \ 1/2 1 r e\ 1/4 1 

The rest follows by re-doing the estimates we did in the case of a dispersion relation of the 
form oj = c s k. Notice that in this case the curvature perturbation ( is of order 

/ TT \ 5/4 -. 

C~#*I M ~(-^J ^~ ir5 H ~ io"V /5 m 2 , (31) 

and the first constraint in (28) becomes 

M 4 > 10 2 ^| . (32) 
Notice also that the cubic operator in M2 generates a three-point function of order 



£2 



H\ 1/A 1,1 

lO- 1 -^, • (33) 



M 2 J 7 7 /2 7 16/5 



Constraints on the three-point function [7, 16] imply 7 > 5. Even though the operator 
generating a four-point function is still 7r 4 the wavefunctions are different in the two cases so 
the induced shape of the four-point function in Fourier space will be a different. In practice 
the difference is not expected to be large, as is the case for the the three-point function [7]. 
On the other hand, it is easy to estimate that the operator 7r 4 that in unitary gauge is 
associated with the operator (5g 00 ) 4 will generate unitary-gauge operators of the form SK 4 
and (5g 00 ) 2 5K 2 , where for simplicity we have suppressed the spacetime indexes which can 
be contracted in all possible ways. These operators will appear with coefficients that make 
them as important as (5g 00 ) 4 . Notice that the scaling dimension for all of these operators is 
the same as the dispersion relation is u oc k 2 . Upon reinsertion of ir, these operators induce 
couplings of the form (<9 2 7r) 4 and 7r 2 (d 2 ii) 2 where the di are spatial derivatives and the indexes 
are contracted in all possible ways. The operators of the form Sg 00 5K 3 and (5g 00 ) 3 SK (with 
spacetime indexes contracted in the several possible ways) are not generate in a relevant way 
because they violate the symmetry n — > —it, t — > —t. On the other hand, we could have 
started directly from these operators of the form 5g 00 SK 3 to find that all the former ones are 
generated with similar size. Summarizing, we see that we can have a total of twelve different 
quartic operators in it involving either one time-derivative or two space derivatives acting 
on each it with various possible tensorial contractions. These twelve different operators have 
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comparable shapes and comparable sizes. This is only possible in the near de Sitter models 
with dispersion relation of the form u = / y 2 k 2 /M 2 . 

In the de Sitter limit and if the dispersion relation is of the form u = c s k (c 2 = HMf / (4M 4 )) 
the unitary gauge operator Sg^SK^^ sets the dispersion relation at horizon crossing. The 
condition for the four-point function induced by the operator (Sg 00 ) A to dominate reads: 

^ S -7% ■ (34) 
H 2 c s ' 

Notice that the three-point function in this case goes as (/c 2 [7], where ( is given by [7] 



M 2 J cl 12 

The constraint in (34) becomes Mf > M^/^c 2 ) in order for the four-point function to be 
the dominant signal. It is consistent to assume that the operators of the form (SK^^) 2 and 
8K U ^K^y that induce a dispersion relation of the form u oc k 2 have a coefficient such that 
they are important only at the unitarity bound induced by the 7r 4 operator. In this case 
the dispersion relation is linear in k up to the cutoff and it is easy to see that the quartic 
operators involving the extrinsic curvature are now generated at a negligible level. In this 
case, the shape of the four-point function is just given by the operator 7r 4 and its shape is 
identical to that in the not-near de Sitter models. If one were to include larger coefficients for 
the operators (5i^%) 2 and 8K V '^SK^^ then one would continuously interpolate to the case of 
a dispersion relation of the form u oc k 2 . 



4 Conclusions 

We have shown that in single-clock inflation it is possible to generate a four-point function 
that could be detected even in the absence of a detectable three-point function. This is 
possible by imposing an approximate continuos shift symmetry and an approximate parity 
symmetry in the Effective Lagrangian for the inflationary fluctuations. For models away from 
the near de Sitter limit, there is a unique interaction operator that gives rise to the large 
four-point function: 7i~ 4 . This results in only one possible shapes for the four-point function. 
Also in the models close to de Sitter the four-point function can be the leading signal, and in 
this case there is more that one possible shape. This motivates us to undertake an analysis 
of the WMAP data in search for such a signal [21]. 
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